3.50 \(\int \frac{d+e x^4}{x^3 \left (a+b x^4+c x^8\right )} \, dx\)

Optimal. Leaf size=199 \[ -\frac{\sqrt{c} \left (\frac{b d-2 a e}{\sqrt{b^2-4 a c}}+d\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} a \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (d-\frac{b d-2 a e}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{2 \sqrt{2} a \sqrt{\sqrt{b^2-4 a c}+b}}-\frac{d}{2 a x^2} \]

[Out]

-d/(2*a*x^2) - (Sqrt[c]*(d + (b*d - 2*a*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sq
rt[c]*x^2)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a*Sqrt[b - Sqrt[b^2 - 4*a*c]
]) - (Sqrt[c]*(d - (b*d - 2*a*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)
/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a*Sqrt[b + Sqrt[b^2 - 4*a*c]])

_______________________________________________________________________________________

Rubi [A]  time = 0.696665, antiderivative size = 199, normalized size of antiderivative = 1., number of steps used = 5, number of rules used = 4, integrand size = 25, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.16 \[ -\frac{\sqrt{c} \left (\frac{b d-2 a e}{\sqrt{b^2-4 a c}}+d\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} a \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (d-\frac{b d-2 a e}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{2 \sqrt{2} a \sqrt{\sqrt{b^2-4 a c}+b}}-\frac{d}{2 a x^2} \]

Antiderivative was successfully verified.

[In]  Int[(d + e*x^4)/(x^3*(a + b*x^4 + c*x^8)),x]

[Out]

-d/(2*a*x^2) - (Sqrt[c]*(d + (b*d - 2*a*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sq
rt[c]*x^2)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a*Sqrt[b - Sqrt[b^2 - 4*a*c]
]) - (Sqrt[c]*(d - (b*d - 2*a*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)
/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a*Sqrt[b + Sqrt[b^2 - 4*a*c]])

_______________________________________________________________________________________

Rubi in Sympy [A]  time = 84.805, size = 201, normalized size = 1.01 \[ - \frac{\sqrt{2} \sqrt{c} \left (2 a e - b d + d \sqrt{- 4 a c + b^{2}}\right ) \operatorname{atan}{\left (\frac{\sqrt{2} \sqrt{c} x^{2}}{\sqrt{b + \sqrt{- 4 a c + b^{2}}}} \right )}}{4 a \sqrt{b + \sqrt{- 4 a c + b^{2}}} \sqrt{- 4 a c + b^{2}}} + \frac{\sqrt{2} \sqrt{c} \left (2 a e - b d - d \sqrt{- 4 a c + b^{2}}\right ) \operatorname{atan}{\left (\frac{\sqrt{2} \sqrt{c} x^{2}}{\sqrt{b - \sqrt{- 4 a c + b^{2}}}} \right )}}{4 a \sqrt{b - \sqrt{- 4 a c + b^{2}}} \sqrt{- 4 a c + b^{2}}} - \frac{d}{2 a x^{2}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  rubi_integrate((e*x**4+d)/x**3/(c*x**8+b*x**4+a),x)

[Out]

-sqrt(2)*sqrt(c)*(2*a*e - b*d + d*sqrt(-4*a*c + b**2))*atan(sqrt(2)*sqrt(c)*x**2
/sqrt(b + sqrt(-4*a*c + b**2)))/(4*a*sqrt(b + sqrt(-4*a*c + b**2))*sqrt(-4*a*c +
 b**2)) + sqrt(2)*sqrt(c)*(2*a*e - b*d - d*sqrt(-4*a*c + b**2))*atan(sqrt(2)*sqr
t(c)*x**2/sqrt(b - sqrt(-4*a*c + b**2)))/(4*a*sqrt(b - sqrt(-4*a*c + b**2))*sqrt
(-4*a*c + b**2)) - d/(2*a*x**2)

_______________________________________________________________________________________

Mathematica [C]  time = 0.0741951, size = 89, normalized size = 0.45 \[ -\frac{\text{RootSum}\left [\text{$\#$1}^8 c+\text{$\#$1}^4 b+a\&,\frac{\text{$\#$1}^4 c d \log (x-\text{$\#$1})-a e \log (x-\text{$\#$1})+b d \log (x-\text{$\#$1})}{2 \text{$\#$1}^6 c+\text{$\#$1}^2 b}\&\right ]}{4 a}-\frac{d}{2 a x^2} \]

Antiderivative was successfully verified.

[In]  Integrate[(d + e*x^4)/(x^3*(a + b*x^4 + c*x^8)),x]

[Out]

-d/(2*a*x^2) - RootSum[a + b*#1^4 + c*#1^8 & , (b*d*Log[x - #1] - a*e*Log[x - #1
] + c*d*Log[x - #1]*#1^4)/(b*#1^2 + 2*c*#1^6) & ]/(4*a)

_______________________________________________________________________________________

Maple [B]  time = 0.025, size = 365, normalized size = 1.8 \[ -{\frac{c\sqrt{2}d}{4\,a}\arctan \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}-{\frac{c\sqrt{2}e}{2}\arctan \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}+{\frac{c\sqrt{2}bd}{4\,a}\arctan \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}+{\frac{c\sqrt{2}d}{4\,a}{\it Artanh} \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}-{\frac{c\sqrt{2}e}{2}{\it Artanh} \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}+{\frac{c\sqrt{2}bd}{4\,a}{\it Artanh} \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}-{\frac{d}{2\,a{x}^{2}}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  int((e*x^4+d)/x^3/(c*x^8+b*x^4+a),x)

[Out]

-1/4*c/a*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x^2*2^(1/2)/((b+(-4*a
*c+b^2)^(1/2))*c)^(1/2))*d-1/2*c/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/
2))*c)^(1/2)*arctan(c*x^2*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*e+1/4*c/a/(-
4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x^2*2^(1/2)/(
(b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b*d+1/4*c/a*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)
^(1/2)*arctanh(c*x^2*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*d-1/2*c/(-4*a*c+
b^2)^(1/2)*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x^2*2^(1/2)/((-b+
(-4*a*c+b^2)^(1/2))*c)^(1/2))*e+1/4*c/a/(-4*a*c+b^2)^(1/2)*2^(1/2)/((-b+(-4*a*c+
b^2)^(1/2))*c)^(1/2)*arctanh(c*x^2*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b*
d-1/2*d/a/x^2

_______________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \[ -\frac{\int \frac{{\left (c d x^{4} + b d - a e\right )} x}{c x^{8} + b x^{4} + a}\,{d x}}{a} - \frac{d}{2 \, a x^{2}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate((e*x^4 + d)/((c*x^8 + b*x^4 + a)*x^3),x, algorithm="maxima")

[Out]

-integrate((c*d*x^4 + b*d - a*e)*x/(c*x^8 + b*x^4 + a), x)/a - 1/2*d/(a*x^2)

_______________________________________________________________________________________

Fricas [A]  time = 0.834555, size = 3742, normalized size = 18.8 \[ \text{result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate((e*x^4 + d)/((c*x^8 + b*x^4 + a)*x^3),x, algorithm="fricas")

[Out]

1/4*(sqrt(1/2)*a*x^2*sqrt(-(a^2*b*e^2 + (b^3 - 3*a*b*c)*d^2 - 2*(a*b^2 - 2*a^2*c
)*d*e + (a^3*b^2 - 4*a^4*c)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^4 - 2*a*b^2*c +
a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)*d^2*e^2)/(a^6*b
^2 - 4*a^7*c)))/(a^3*b^2 - 4*a^4*c))*log((3*a*b^2*c*d^2*e^2 - 3*a^2*b*c*d*e^3 +
a^3*c*e^4 + (b^2*c^2 - a*c^3)*d^4 - (b^3*c + a*b*c^2)*d^3*e)*x^2 + 1/2*sqrt(1/2)
*((b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*d^3 - (3*a*b^4 - 13*a^2*b^2*c + 4*a^3*c^2)*d^2
*e + 3*(a^2*b^3 - 4*a^3*b*c)*d*e^2 - (a^3*b^2 - 4*a^4*c)*e^3 - ((a^3*b^4 - 6*a^4
*b^2*c + 8*a^5*c^2)*d - (a^4*b^3 - 4*a^5*b*c)*e)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4
- (b^4 - 2*a*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a
^3*c)*d^2*e^2)/(a^6*b^2 - 4*a^7*c)))*sqrt(-(a^2*b*e^2 + (b^3 - 3*a*b*c)*d^2 - 2*
(a*b^2 - 2*a^2*c)*d*e + (a^3*b^2 - 4*a^4*c)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^
4 - 2*a*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)
*d^2*e^2)/(a^6*b^2 - 4*a^7*c)))/(a^3*b^2 - 4*a^4*c))) - sqrt(1/2)*a*x^2*sqrt(-(a
^2*b*e^2 + (b^3 - 3*a*b*c)*d^2 - 2*(a*b^2 - 2*a^2*c)*d*e + (a^3*b^2 - 4*a^4*c)*s
qrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^4 - 2*a*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2
*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)*d^2*e^2)/(a^6*b^2 - 4*a^7*c)))/(a^3*b^2 - 4*
a^4*c))*log((3*a*b^2*c*d^2*e^2 - 3*a^2*b*c*d*e^3 + a^3*c*e^4 + (b^2*c^2 - a*c^3)
*d^4 - (b^3*c + a*b*c^2)*d^3*e)*x^2 - 1/2*sqrt(1/2)*((b^5 - 5*a*b^3*c + 4*a^2*b*
c^2)*d^3 - (3*a*b^4 - 13*a^2*b^2*c + 4*a^3*c^2)*d^2*e + 3*(a^2*b^3 - 4*a^3*b*c)*
d*e^2 - (a^3*b^2 - 4*a^4*c)*e^3 - ((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*d - (a^4*
b^3 - 4*a^5*b*c)*e)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^4 - 2*a*b^2*c + a^2*c^2)
*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)*d^2*e^2)/(a^6*b^2 - 4*a
^7*c)))*sqrt(-(a^2*b*e^2 + (b^3 - 3*a*b*c)*d^2 - 2*(a*b^2 - 2*a^2*c)*d*e + (a^3*
b^2 - 4*a^4*c)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^4 - 2*a*b^2*c + a^2*c^2)*d^4
+ 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)*d^2*e^2)/(a^6*b^2 - 4*a^7*c)
))/(a^3*b^2 - 4*a^4*c))) + sqrt(1/2)*a*x^2*sqrt(-(a^2*b*e^2 + (b^3 - 3*a*b*c)*d^
2 - 2*(a*b^2 - 2*a^2*c)*d*e - (a^3*b^2 - 4*a^4*c)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4
 - (b^4 - 2*a*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 -
a^3*c)*d^2*e^2)/(a^6*b^2 - 4*a^7*c)))/(a^3*b^2 - 4*a^4*c))*log((3*a*b^2*c*d^2*e^
2 - 3*a^2*b*c*d*e^3 + a^3*c*e^4 + (b^2*c^2 - a*c^3)*d^4 - (b^3*c + a*b*c^2)*d^3*
e)*x^2 + 1/2*sqrt(1/2)*((b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*d^3 - (3*a*b^4 - 13*a^2*
b^2*c + 4*a^3*c^2)*d^2*e + 3*(a^2*b^3 - 4*a^3*b*c)*d*e^2 - (a^3*b^2 - 4*a^4*c)*e
^3 + ((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*d - (a^4*b^3 - 4*a^5*b*c)*e)*sqrt(-(4*
a^3*b*d*e^3 - a^4*e^4 - (b^4 - 2*a*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^
3*e - 2*(3*a^2*b^2 - a^3*c)*d^2*e^2)/(a^6*b^2 - 4*a^7*c)))*sqrt(-(a^2*b*e^2 + (b
^3 - 3*a*b*c)*d^2 - 2*(a*b^2 - 2*a^2*c)*d*e - (a^3*b^2 - 4*a^4*c)*sqrt(-(4*a^3*b
*d*e^3 - a^4*e^4 - (b^4 - 2*a*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e -
 2*(3*a^2*b^2 - a^3*c)*d^2*e^2)/(a^6*b^2 - 4*a^7*c)))/(a^3*b^2 - 4*a^4*c))) - sq
rt(1/2)*a*x^2*sqrt(-(a^2*b*e^2 + (b^3 - 3*a*b*c)*d^2 - 2*(a*b^2 - 2*a^2*c)*d*e -
 (a^3*b^2 - 4*a^4*c)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^4 - 2*a*b^2*c + a^2*c^2
)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)*d^2*e^2)/(a^6*b^2 - 4*
a^7*c)))/(a^3*b^2 - 4*a^4*c))*log((3*a*b^2*c*d^2*e^2 - 3*a^2*b*c*d*e^3 + a^3*c*e
^4 + (b^2*c^2 - a*c^3)*d^4 - (b^3*c + a*b*c^2)*d^3*e)*x^2 - 1/2*sqrt(1/2)*((b^5
- 5*a*b^3*c + 4*a^2*b*c^2)*d^3 - (3*a*b^4 - 13*a^2*b^2*c + 4*a^3*c^2)*d^2*e + 3*
(a^2*b^3 - 4*a^3*b*c)*d*e^2 - (a^3*b^2 - 4*a^4*c)*e^3 + ((a^3*b^4 - 6*a^4*b^2*c
+ 8*a^5*c^2)*d - (a^4*b^3 - 4*a^5*b*c)*e)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^4
- 2*a*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)*d
^2*e^2)/(a^6*b^2 - 4*a^7*c)))*sqrt(-(a^2*b*e^2 + (b^3 - 3*a*b*c)*d^2 - 2*(a*b^2
- 2*a^2*c)*d*e - (a^3*b^2 - 4*a^4*c)*sqrt(-(4*a^3*b*d*e^3 - a^4*e^4 - (b^4 - 2*a
*b^2*c + a^2*c^2)*d^4 + 4*(a*b^3 - a^2*b*c)*d^3*e - 2*(3*a^2*b^2 - a^3*c)*d^2*e^
2)/(a^6*b^2 - 4*a^7*c)))/(a^3*b^2 - 4*a^4*c))) - 2*d)/(a*x^2)

_______________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \[ \text{Timed out} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate((e*x**4+d)/x**3/(c*x**8+b*x**4+a),x)

[Out]

Timed out

_______________________________________________________________________________________

GIAC/XCAS [A]  time = 0.495389, size = 1, normalized size = 0.01 \[ \mathit{Done} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate((e*x^4 + d)/((c*x^8 + b*x^4 + a)*x^3),x, algorithm="giac")

[Out]

Done